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Abstract 

Original Whitham’s method of derivation of modulation equations is applied to systems whose dynamics is 
described by a perturbed Korteweg-de Vries equation. Two situations are distinguished: (i) the perturbation 
leads to appearance of right-hand sides in the modulation equations so that they become non-uniform; (ii) 
the perturbation leads to modification of the matrix of Whitham velocities. General form of Whitham 
modulation equations is obtained for each case. The essential difference between them is illustrated by an 
example of so-called ‘generalized Korteweg-de Vries equation’. Method of finding steady-state solutions of 
perturbed Whitham equations in the case of dissipative perturbations is considered. 
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Dedicated to fiftieth anniversary of publication of 
Whitham’s ‘the Paper’. 

1. Introduction 

In his seminal paper [i| Whitham introduced into 
nonlinear wave theory several fundamental ideas 
which formed the basis for development of a vast 
theory called now Whitham theory. First, he gen¬ 
eralized the idea of slow evolution of envelopes of 
linear harmonic wave trains (‘wave packets’) to de¬ 
scription of evolution of nonlinear modulated wave 
trains whose dynamics is governed by nonlinear 
wave equations. This idea implies that in the 
problem under consideration there are two different 
scales of space and time: the field variables u{x,t) 
of the nonlinear ‘carrier’ wave oscillate at the scales 
of wavelength L and period T, whereas such param¬ 
eters of the wave as, e.g., wavelength L, amplitude 
a, phase velocity V, etc., change slowly at the space 
scale a; » L and the time scale t^ T. This leads 
to the second idea of averaging of the conservation 
laws of the evolution equation over fast local oscilla¬ 
tions analogously to the Krylov-Bogoliubov averag¬ 
ing technique developed in the theory of nonlinear 
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vibrations. However, in contrast with dynamical 
time-dependent systems, now the field variables de¬ 
pend on time and one (or more) space coordinates 
and, as a result of averaging of conservation laws, 
Whitham obtained the system of first order par¬ 
tial differential equations now called Whitham equa¬ 
tions. Whitham compared this approach with tran¬ 
sition from ‘microscopic’ description of gas dynam¬ 
ics to averaged hydrodynamic description (“Indeed, 
the present work is in much the same spirit as the 
derivation of continuum fluid mechanics from ki¬ 
netic theory.”) and this suggested the third idea of 
application of the averaged equations to description 
of such physical phenomena as, for example, water 
undular .bores and collisionless shocks in plasma. 
At last, as the forth idea, Whitham supposed that 
his modulation equations can be transformed, by 
analogy with compressible fluid dynamics, to the 
diagonal Riemann form and he realized this idea 
by means of very skillful calculations for the case 
of modulated nonlinear ‘cnoidal’ wave whose evolu¬ 
tion is governed by the Korteweg-de Vries ^dV) 
equation. Thus, Whitham formulated in [ij the 
general method for studying modulated nonlinear 
waves and illustrated fruitfulness of his approach 
by important nontrivial examples. Richness of ideas 
introduced in [l[ has been spectacularly confirmed 
by impressive development of the Whitham theory 
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in the past 50 years. 

The first important contribution into the 
Whitham theory after appearance of the paper 
[H was done by Gurevich and Pitaevskii [J who 
showed that a collisionless shock (now commonly 
called dispersive shock wave (DSW)) described by 
the KdV equation can be represented as an expand¬ 
ing oscillating structure which can be approximated 
by a modulated cnoidal wave whose evolution is 
governed by the Whitham equations. At one its 
edge the dispersive shock approaches to a soliton 
train and at the opposite edge it tends to a small 
amplitude harmonic wave. Gurevich and Pitaevskii 
studied self-similar solutions of the Whitham equa¬ 
tions for a typical examples of evolutions of an ini¬ 
tial step-like distribution and of the general wave 
breaking situation. 

Analytical theory developed by Gurevich and 
Pitaevskii was based on a specific diagonal form 
of modulation equations obtained by Whitham for 
modulated KdV cnoidal waves. However, such a 
form for other nonlinear wave equations was not 
known and it was not easy to find it by the di¬ 
rect method used by Whitham. Actually, as it be¬ 
came clear later, Whitham had succeeded in find¬ 
ing the Riemann invariants for the KdV equation 
case because this equation belongs to a very special 
class of so-called ‘completely integrable equations’ 
whose solutions can be found by the inverse scat¬ 
tering transform (1ST) method discovered indepen¬ 
dently of the Whitham theory [ 1 , 0 , 0 • It turned 
out that this method generalized on quasi-periodic 
situations 00 yields quasi-periodic solutions of 
the KdV equation which are parameterized directly 
by Riemann invariants having in this case very sim¬ 
ple mathematical meaning: they are the edge points 
of gaps in the spectrum of the linear (Schrodinger) 
equation related with the KdV equation in the 1ST 
method. As a result, multi-phase averaging method 
for the quasi-periodic solutions of the KdV equation 
was developed by Flaschka, Forest and McLaughlin 
0 where the Whitham equations were derived in 
Riemann diagonal form for 2A^ -|- 1 dependent vari¬ 
ables (Riemann invariants), N being the number of 
phases (TV = 1 for the sirnplest cnoidal wave case 
considered by Whitham in 0). The general method 
of derivation of the Whitham equations for a wide 
class of completely integrable equations was sug¬ 
gested by Krichever 0 . The generalized hodograph 
method of integration of diagonal Whitham equa¬ 
tions was developed by Tsarev [ 10 . This progress 
in mathematical theory of integrability of nonlin¬ 


ear wave equations and of corresponding Whitham 
modulation equations has led to a number of appli¬ 
cations to physical problems related with formation 
of DSWs and deeper understanding of qualitative 
properties of this phenomenon. At the same time, 
it became quite desirable to extend the theory of 
DSWs on situations often met in physical appli¬ 
cations when wave motion is not described by the 
completely integrable equations. 

From physical point of view, it seems clear that 
the phenomenon of formation of DSWs is related 
with effects of dispersion in nonlinear wave systems 
and is not conditioned by complete integrability of 
corresponding evolution equations. Actually, the 
Whitham theory was developed in [ij in very gen¬ 
eral setting under supposition of existence of peri¬ 
odic solutions of evolution equations and only appli¬ 
cation of this theory to the KdV cnoidal wave was 
related implicitly in this paper with the complete 
integrability of the KdV equation. Therefore for 
treatment of DSWs in general situation one should 
resort to analysis of Whitham equations in a non¬ 
diagonal form when they do not have Riemann in¬ 
variants and cannot be integrated by Tsarev’s gen¬ 
eralized hodograph method. Such an analysis was 
done by El 11| in an important particular case of 
evolution of an initial step-like pulse whose dynam¬ 
ics is governed by non-dissipative nonlinear wave 
equations. This ingenious method has found a num¬ 
ber of interesting applications in which the problem 
can be reduced to the study of evolution of step-like 
pulses. 

Another typical situation appears when the evo¬ 
lution equation differs little from an integrable one. 
For example, such a difference can appear due to 
small dissipation effects or weak non-uniformity of 
the medium through which the wave prcmagates. 
As was indicated already by Whitham in [l|, these 
effects lead to modulation of nonlinear periodic 
waves and can be considered in framework of the 
averaged modulation equations. In such situations, 
the Whitham modulation equations can be modi- 
hed by perturbations in two possible ways: (i) the 
equations for Riemann invariants Xi{x,t) of unper¬ 
turbed equations have now ‘right-hand sides’ de¬ 
pending on Ai, that is these equations become non- 
uniform; (ii) the additional terms caused by per¬ 
turbations contribute to fluxes of the conserved 
quantities leading to appearance in the Whitham 
equations of terms proportional to the derivatives 
dXj/dx, that is the matrix of ‘velocities’ becomes 
non-diagonal and Xi are not Riemann invariants 
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anymore, but the non-diagonal terms as well as cor¬ 
rections to the diagonal ones are small and the mod¬ 
ulation equations remain uniform. Of course, one 
can imagine situations when both types of correc¬ 
tions appear in perturbed Whitham equations. 

So far, mainly the first type of corrections has 
been considered. Physically, such corrections ap¬ 
pear very naturally when irreversible processes 
are taken into account. The Whitham averag¬ 
ing method for systems with small dissipation was 


developed by Jimenez and Whitham 12| in gen¬ 


eral form without transition to Riemann invari¬ 
ants. Whitham equations for fV-phase KdV wave 
trains in presence of small perturbations were de¬ 
rived in [l3|, however in a form not convenient 
enough for applications. More practical and in¬ 
structive example of one-phase modulated KdV 
wave trains with account of small Burgers viscos¬ 
ity was considered by Gurevich and Pitaevskii 0 
and by Avilov, Krichever and Novikov [l^ (earlier 
the steady-state solution of this problem had been 
studied by Johnson [16| by a direct perturbation 
technique). They derived the Whitham equations 
for the Riemann invariants Xi of unperturbed prob¬ 
lem and showed that small Burgers viscosity results 
in non-zero right-hand sides of Whitham equations 
which provide additional contribution into evolu¬ 
tion of these modulation parameters Xi. The anal¬ 
ysis presented there showed that although the per¬ 
turbation is small compared with the main terms in 
the KdV equation, this does not mean that its con¬ 
tribution into evolution of the Riemann invariants 
is also small. Indeed, in this case the perturbation 
should be compared with a small parameter which 
characterizes ‘slowness’ of modulation rather than 
with terms which determine fast oscillations of the 
cnoidal wave. If the cnoidal wave is not modulated 
at all, then dissipative terms make the only contri¬ 
bution into changes of the Riemann invariants and 
determine slow evolution of a uniform cnoidal wave 
(see, e.g., [ 13 , 0 ). Effects of non-local damping 
were considered by Gurevich and Pitaevskii in [19l | 
and more general forms of local dissipation were 
considered by Myint and Grimshaw in [0. Quite 
general approach applicable to the Ablowitz-Kaup- 
Newell-Segur (AKNS) class 21| of completely in- 
tegrable equations was developed by the author 
for non-perturbed and perturbed [^ cases. 
In combination with simplified version |24| of the 
finite-gap integration method which yields the pe¬ 
riodic solutions in a ‘real’ form not-constrained by 
any additional ‘reality conditions’, this approach 


turned out to be quite effective and it has found sev¬ 
eral non-trivial applications including propagation 
of KdV wave trains through a non-uniform medium 
(see, e.g., 


However, the approach described above is not 
applicable to situations when perturbations change 
the matrix of Whitham velocities although such sit¬ 
uations are encountered quite often. For example, 
if in the perturbed KdV equation 


Ut + QuUx + Uxxx = -RM (1) 

the perturbation term has the form R[u\ = 
eF'{u)ux,(€ <C 1), then the general formulae ob¬ 
tained in [i^, 0 lead to vanishing right-hand sides 
in the ‘perturbed’ Whitham equations what means 
that such a perturbation belongs to the type (ii) and 
a different perturbation scheme should be devel¬ 
oped for finding the corrected matrix of Whitham 
velocities. This difference between two types of per¬ 
turbations is clearly illustrated by a simple example 
of perturbation i?[u] = const-u^Ua, when ([T]) reduces 
to the so-called Gardner equation. This equation is 
also completely integrable, the corresponding Rie¬ 
mann invariants and Whitham equations can be ob¬ 
tained without any approximations (see 00), 
and they do not reduce to appearance of the right- 
hand sides in Whitham equations in the KdV limit 
of the Gardner equation. 

Thus, we arrive at the problem of derivation of 
the approximate Whitham equations for two differ¬ 
ent situations when either the Whitham equations 
acquire the right-hand side terms, or the Whitham 
velocities are modified by perturbations. Here we 
shall confine ourselves to a simple example of the 
KdV equation CD under supposition that the per¬ 
turbation term is small, 

|i?[u]| <C min{u^/L, |u|/L^}. (2) 

Hence, we can approximate locally the solution of 
CD by the cnoidal wave solution of unperturbed 
KdV equation and apply the original method of 
Whitham [i| to this more general situation. In the 
next section we shall illustrate the method by its 
application to the already studied earlier situation 
of perturbations of type (i) and then generalize it 
to perturbations of type (ii). In section 3 we shall 
show that a specific structure of perturbation terms 
leads to a simple method of finding the steady-state 
solutions of the Whitham equations. We conclude 
by the remark that the direct Whitham approach 
to obtaining the modulation equations is effective 
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enough and it can be successfully used for studying 
quite complicated nonlinear wave problems. 


2. Perturbed Whitham equations 

Traveling wave solution of the unperturbed KdV 
equation 

Ut + &UUx + ^XXX - 0 ( 3 ) 

is obtained by means of a simple substitution u = 
m(^), ^ = x — Vt, so that after two obvious integra¬ 
tions we get the ordinary differential equation 

\ul = f{u), (4) 

where f{u) is a third degree polynomial (a > /3 > 
7). 

f{u) = -A + Bu+^Vu^-u^ 

= -(w-a)(M-/3)(u-7), 

where V, and B are the integration constants re¬ 
lated with the zeros a, /?, 7 of the polynomial f{u) 
by the formulae 

V = 2{a + l3 + j), A = -al3'y, 

B = —{af3 + + 7 a). 

In a standard way the solution of Eq. o can be 
expressed in terms of Jacobi elliptic sinus function 

u{x,t) = a- {a- /3)sn2(^/ {a - j)/2^,m), (7) 


where the parameter m is equal to 


a — 7 


( 8 ) 


This is a periodic solution of the KdV equation 
and its wavelength is given by the formula 


1 


du 


L = -= d£,= (h — 


1 ^ du 

\f 2 J a//(m) 


= 2 


(9) 


a — 7 


■K{m), 


where K{m) is the complete elliptic integral of the 
first kind. We have introduced in m the wavenum¬ 
ber k = 1/L and hence at a given point x the wave 
oscillates with the frequency uj = kV, so that the 
solution © depends on the phase 9 = kx — tot. 

In a modulated wave the parameters V, A, 
or, equivalently, a, j5, 7 , become slow functions of 
X and t, that is they change little in one wavelength 


L and one period T ~ l/w. However, according to 
Whitham, the solution of the (perturbed or unper¬ 
turbed) KdV equation can be approximated locally 
by the expression © where a{x,t), ld{x,t), 
are considered now as slow functions of x and t and 
the phase 9 = kx — ujt is replaced now by a general 
dependence 9{x,t)- Then the wavenumber and the 
frequency are defined as 

k = 9x, UJ = -9t, (10) 

and, hence, they must satisfy the compatibility con¬ 
dition kt + uJx =0 which has the meaning of con¬ 
servation of ‘number of waves’ [ij. In a slowly mod¬ 
ulated wave both k and w = kV are expressed in 
terms of the slow parameters a, /3, 7 (see Eqs. (O 
and (O) so that we arrive at the equation for these 
parameters, 

kt + {kV)x = 0. (11) 

All that is applied to any modulated KdV wave 
train and the modulation can be caused either by 
a non-uniform initial condition of by a perturba¬ 
tion term in ©• Whitham discussed in [l| the first 
situation only and we wish here to generalize his ap¬ 
proach to perturbed KdV equations. As was indi¬ 
cated in Introduction, we have to distinguish in this 
case two different situations which, as we shall see, 
can be formulated more precisely as follows: (i) nei¬ 
ther R nor uR are space derivatives, (ii) R and/or 
uR can be represented as space derivatives of other 
functions (say, R = Qi^x and/or uR = Q 2 ,x)- We 
shall call the first situation as a non-gradient per¬ 
turbation and the second one as a gradient pertur¬ 
bation, and we shall begin with discussion of the 
non-gradient perturbations. 


2.1. Whitham eguations for the case of non¬ 
gradient perturbations 

In addition to CD, we need two more equa¬ 
tions for three parameters a{x,t), j3{x,t), j{x,t) 
and, following Whitham [l[, we assume that they 
can be obtained by means of averaging the conser¬ 
vation laws of the KdV equation o, 

“t” (3w J- Uxx)x — R, (12) 

+ (2m^ -I- UUxx - ^ul)x = uR. 

Averaging is defined as taking a mean value of an 
expression V along the wavelength. 


1 


Vdx = k 


^du k 
V — = —1= 
Ux V2 


Vdu 


( 13 ) 
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where integration is taken over the whole cycle of 
oscillation of u. Whitham averaged densities and 
fluxes of the conservation laws (1121) (with R = 0) 
and obtained two additional equations for the slow 
variables. If we average the right-hand sides of 
Eqs. (fT^ according to the rule (fT^ then we obtain 
the perturbed Whitham equations with the right- 
hand sides. The averaged conservation laws m 
take the form 


and substitute (fTCl) into (ITT|) to obtain the complete 
set of the Whitham modulation equations: 


DWa 

Dt 

DWb 

Dt 

DWv 

Dt 


Wa-^, 

ox 

Wa-^-Wa{R), 

dA 

Wa^-Wa{uR). 


(19) 


{u)t + (3u^ + Uxx)x — {R), 

{W)t + (2m^ -f uuxx - hul)x = (uR) 


(14) 


and the condition that R and uR are not x- 
derivatives yields, generally speaking, non-zero 
right-hand sides in these equations. The derivatives 
u^, Uxx can be excluded with the use of Eqs. o, 
(O: = 2/(m), Uxx = f'{u) = B + Vu — 3m^, and, 

as Whitham indicated, it is convenient to represent 
the averaged expressions in terms of a single ‘action 
function’ 


W{A, B,V) = -V2ji du 


= —V 2 f \ —A + Bu \Vu^ — du, 


(15) 


and its derivatives with respect to A, B, V. Indeed, 
the wavelength dOj can be written as 


L = i = ±l 


dW ... 

= j^^Wa (16) 


k 72/ /7()I) dA 
and the necessary averages are given by 

/I 2 \ k [{u^/2)du 
{^u^) = —;= <p ! A, , = -kWv- 


V2J yW) 


(17) 


As a result we obtain the averaged conservation 
laws 


{-kWB)t + {B- kVWB)x = (R), 
{-kWv)t + {A- kVWv)x = (uR), 


which can be simplified with the use of Eq. (HH). 
Besides that, we introduce the ‘long derivative’ 


_D _ d 

Dt dt dx 


Naturally, they differ from the original Whitham 
equations [i| only by the terms with (R) and {uR). 

A remarkable discovery of Whitham was that the 
unperturbed modulation equations can be trans¬ 
formed “after considerable manipulation” to the di¬ 
agonal (Riemann) form and the Riemann invariants 
are expressed in terms of zeros a, / 3 , 7 . In our case 
the same transformation leads again to the diago¬ 
nal form of Whitham equations, however now with 
small right-hand sides, that is the Whitham equa¬ 
tions become non-uniform. The necessary ‘consid¬ 


m 


erable manipulation’ is described in detail in 
and we indicate here briefly the main steps only. 

First, we transform Eqs. m from the variables 
A, B, V to the variables a, j3 ,7 with the use of re¬ 
lationships dH), 


Wa. 


Da 


Dt 
= 2WA(a, 
Da 


Wa,p 


£1 

Dt 

klx), 

Dfd 


Wa. 


D'y 

In 




D'y 


X [(/3 + 7)aa; + (a -I- ^)fix + (a + P)lx] 

-Wa{R), 

X [^7 ■ ax + a^ ■ I3x + al3 ■ jx] - Wa{uR), 


( 20 ) 


where 


WA,a. 


WB,a 


Wv,a. 


1 / du 

Vs f {u- a)Vf{u)' 

1 / udu 

VSJ {u- a)Vf{u)' 

1 {V/2)du 

Vsf {u- a) Vf{uy 


( 21 ) 


and similar expressions can be written for deriva¬ 
tives with respect to /3 and 7 . Next, we multiply the 
first equation (pnj) by p = X(3 + a'y — /Sy, the second 
equation by g = 2a, the third equation by r = 2 


(18) 
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(28) 


and add them; then with the use of the identities 


and the wavelength is given by 


+ If A, = = 0, 

pWA,a + qWB,a + rWv,a 




we obtain 

+ 7) d{f3 + 7 ) 


Dt WA,a dx 
1 Wa 

(a - /3){a - 7 ) Wa.q 


( 22 ) 


(a(i?) — (uR)) 


and similar equation can be obtained for the vari¬ 
ables a + P and a -I- 7 by means of cyclic trans¬ 
position of the parameters a, /3, 7 . At last, we in¬ 
troduce the Riemann invariants of the unperturbed 
KdV equation, 


(23) 


Ai =-i(a-b/3), X2 = + 

A 3 = + 7 ), Ai < A 2 < A 3 , 

and with account of 

Wa 2Wa 2 L 
“ Wa,a3 “ dL/dX3 

and similar formulae for Wa/WA, p and Wa!Wa,' i 
we arrive at the Whitham equation in the form 

dXj ( 0 ) dXj L {{2Xi - Si - u)R) 
dt dx dhjdXi ‘^Y{j^i{Xi - Xj) (24) 

* = 1,2,3, 

where velocities are given by 

and Si = A 1 -I-A 2 -I-A 3 . Here all the variables 
should be parameterized by the Riemann invariants 
Ai, A 2 , A 3 . In particular, the periodic solution of 
the KdV equation takes the form 


l{x, t) — X 3 — X 2 — Ai 


- 2 (A 3 - A 2 ) sn^(-\/A 3 - Ai (x - Vt), m), 
where 

V = —2si = —2(Ai -b A 2 -b A 3 ), 

A 3 — A 2 

" A 3 - A, ’ 


(26) 


(27) 


^ _ 1 _ 2K{m) 
k a/As — Ai 

Substitution of (1^ into (ES|) gives expressions for 
the Whitham velocities in their original form [I|, 

„(o) _ „ 4(A3-Ai)(l-m)A:(TO) 

V-^ — ZSi “T - 


= — 2 si — 

^ E(rn) — {1 — m)K{m) 

4 °) = - 2 si ' ■ ^ 2 )K(m) 


E{m) 

4 (A 3 — A 2)(1 — m)K{m) 


, (29) 


E{m) — K{m) ’ 


where E(rn) is the complete elliptic integral of the 
second kind. For averaging the perturbation terms, 
it is convenient to introduce the variable fj, = (u + 
si )/2 which changes within the interval A 2 < /r < 
A 3 and satisfies the equation 


Ha; = 2 ^J-P{h), 


(30) 


where 


Pip) = ip->i2)ip-X3) 

,3 „ .,2 . 


= H - Sifi +S2P-S3, 


consequently 


1 




<«> = ! 




Qdn 

'a. 


(32) 


These formulae easily reproduce the results found 
in Refs. l3. 15. IT^ with R, = £Uxx corresponding to 
the Burgers viscosity and permitted one to derive 
Whitham equations for the case of i? = F{t)u — 
Git)v? corresponding to shallow water waves over a 
gradual slope with account of bottom friction . 
However, as was mentioned in Introduction, this 
approach fails if averages of R and/or uR vanish 
what happens when these expressions are the x- 
derivatives of some other expressions and, hence, 
they contribute into the fluxes of the conservation 
laws (HU. We shall consider such a situation in the 
next subsection. 

2.2. Whitham equations for the case of gradient 
perturbations 

Here we shall assume that both R and uR are 
space derivatives 

R — Ql,x; uR — 0,2,X. 
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(33) 

























and, hence, they make additional contribution into 
the fluxes in the conservation laws (fT^ . Then all 
divergence terms in m should be treated on the 
same footing and averaging of the conservation laws 
yields 


{3u Uxx'jx — (Ql)a 


+ {2u^ + UUxx - 


2 


)x — {Q2 )x- 


(34) 


Transformations similar to those which were done 
above give instead of Eq. (|TO)) the equations 


Dt ^ dx ^ dx '' 

Dt ^dx ^ dx ' 


(35) 


or, after transition to the variables a, /3, 7 , the 
equations 

Da DB £>7 


— 2Wa(0;x + jBx + 7a:), 

Da DB 


D^ 


X [(/3 + l)(^x + (a + 'y)/3x + (a + /3)jx] 
— Wa{{Qi) aOtx + (Qi)/3/3x + (Ql)77a:)! 
Da 


(36) 


Wv,c 


Dt 


DB £>7 

Wv,P^ + Wv,,^ 


= -WaIPj ■ ax + aj ■ Px + aP ■ 7a;] 

~ Wa{{Q 2)aOtx + {Q2)pfix + {Q2)~f'yx)- 

Their linear combination used above for transition 
to Eq. (1^ now gives 


£)(/3 + 7 ) , Wa d{jB + 7 ) 


Dt 


Wa,, 


dx 

Wa 


{(a(Ql}a - (Q2}a)ax 


{a - /3){a - 7 ) 1Ta,c 

+ (/3(Qi)/3 — {Q2)p)I3x + ( 7 ( 21)7 — {Q2)'f)"yx}- 

(37) 

At last, transformation to the variables Ai, A 2 , A 3 
(see (USD) yields the Whitham equations in the form 


3Ai ( 0 ) dXi 


1 


E ( 2 A - si) 


4n,^a(A,-A,)5£/5A 
d{Qi) 5(22) \ dXk (38) 


dXk 


dXk 


dx ’ 


i = 1,2,3. 


Here it was supposed that both R and uR are gra¬ 
dients of ‘fluxes’ 2i and Q 2 , respectively. If only 
one of the variables R and uR can be represented 
as a gradient, then only one corresponding term 
(2fc) is included in (1551) and the other term must 
be treated as a non-gradient one resulting in the 
right-hand sides of the Whitham equations, as was 
shown in the preceding subsection. 

Generally speaking, the variables Ai, A 2 , A 3 in 
Eqs. (1381) are not Riemann invariants anymore. 
Indeed, we have got a non-diagonal matrix of 
Whitham velocities 


dt 


3 


El 

' dx 


= 0 , 


,( 0 ); 


^13 ’ 


(39) 


where correspond to the perturbation terms. 
They are much smaller than the contributions 
v^^'^Sij into the diagonal ones and we can find the 
characteristic velocities as well as the correspond¬ 
ing eigenvectors in the way similar to the stationary 
perturbation theory well known in quantum me¬ 
chanics (see, e.g., section 38 in [^) as long as 

> max|u-j^|. (40) 

In particular, the characteristic velocities are given 
in our approximation by the formulae 


~ ( 0 ) , 

Vi = vy+ V, 


„( 1 ) 


(41) 


The condition (HU)) indicates that Eqs. (1551) cannot 
be applied to situations with degeneration of two 
Whitham velocities, that is, for example, to the 
edge points of dispersive shock waves. This con¬ 
clusion is confirmed by the following remark. As 
is known (see, e.g., i) , mean value of the variable 
u vanishes at the soliton edge of dispersive shock 
wave according to the law 


(u) oc 


1 

ln(A2 — Ai) 


(42) 


Hence, this mean value has infinite derivatives with 
respect to Riemann invariants Ai and A 2 , and mean 
values (u") have similar singularities at the soliton 
edge. Since typically the perturbation terms are 
expressed in terms of such mean values and their 
derivatives, then the characteristic velocities m 
are also singular at the soliton edge what prevents 
application of Eqs. (1551) to the theory of DSWs. A 
possible method of avoiding this difficulty is dis¬ 
cussed in the following subsection. 
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2.3. Elimination of gradient perturbations from 
Whitham equations 

There are situations when the gradient terms can 
be eliminated from the Whitham equations and the 
difficulty indicated in the preceding subsection can 
be avoided by means of the method used first by 
Marchant and Smyth in application of the ex¬ 
tended KdV (or Gardner) equation 

ut + 6uux + Uxxx = Qeu^Ux (43) 

to formation of undular bores in the resonant flow 
of a fluid over topography. It is supposed that e is a 
small parameter, e <C 1, and Marchant and Smyth 
showed that with accuracy 0 (e) the equation can 
be reduced to the KdV equation 

Ut + WUx + Uxxx = 0 (44) 

by means of a simple substitution 

u = U + e(t7^ -|- Uxx/‘2‘) 
or (45) 

U = u- e{u^ +Uxxl‘i)- 

In fact, the equation (H5)) is completely integrable 
and the Whitham equations can be derived [i^, 
beyond the perturbation theory for any value (and 
sign) of e. Nevertheless, as we shall show, the 
method of substitutions similar to that of Marchant 
and Smyth turns out to be quite useful in a situa¬ 
tion of the ‘generalized KdV equation’ 

Ut + 6uux + Uxxx = 6eF{u)ux, (46) 

where F{u) is a regular function in the region of 
variations of u. 

Let us look at several examples of simple sub¬ 
stitutions and the results of the corresponding 
transformations. Everywhere we neglect the terms 
smaller than the order 0{e). 

A substitution 

U = u - euxx or u = U + eUxx (47) 

leads to an approximate (with accepted here accu¬ 
racy) identity 

Ut UUUx ~\~ Uxxx — 4“ QuUx Uxxx T ^^.eUxUxx- 

Hence, the perturbed KdV equation 

Ut + QuUx + Uxxx =-^‘^^UxUxx (48) 


is reduced to Eq. (HU) by means of the substitution 

(071). 

A substitution 

U = u — eu^ or u = U + eU^ (49) 

transforms in a similar way the perturbed KdV 
equation 

Ut + 6uUx + Uxxx = 6e{u‘^Ux + UxUxx) (50) 

to Eq. {|44j) . Composition of substitutions 011) and 
(l4^ yields the substitution (1451) of Marchant and 
Smyth. In this particular case the generalized KdV 
equation (H51) is reduced to the non-perturbed KdV 
equation (H41) , however such a reduction is generally 
impossible. For example, a substitution 


U = u — eu^ or u = U + eU^ 

(51) 

reduces the equation 


Ut + 6uux + Uxxx = Geu^Ux 

(52) 

to 


Ut + 6UUx + Uxxx = 6e(C/3 + SUUxUxx), 

(53) 

and here the terms in the right-hand side cannot 
be eliminated by additional substitutions. In the 
general case a substitution 

U = u — eF{u) or u = U + eF{U) 

(54) 

transforms (H51) to 


Ut + GUUx + Uxxx 

= &e{F'"(U)Ul + SF"iU)UxUxx)- 

(55) 


Thus, these substitutions transform a perturbed 
equation to another perturbed equation, however 
there is an important difference between the initial 
and reduced forms: the perturbation in Eq. 051) is 
a gradient one whereas in Eq. ()55l) it is non-gradient 
and therefore these perturbations should be treated 
by different types of Whitham equations discussed 
above. 

Another important feature of the perturbation 
terms in (155|) is that they do not contribute into 
the right-hand sides of the Whitham equations (1241) 
which coincide, hence, with unperturbed Whitham 
equations for the KdV equation. Let us demon¬ 
strate this for the averaged value (i?). Simple inte¬ 
gration by parts shows that 

j F"UxUxxdx = -\j 


F"'{U)Uldx 


that is contribution of the second term in the right- 
hand side of Eq. ((551) reduces after averaging (up to 
factor—1/2) to the contribution of the first term. 
But its average vanishes, as shows a simple calcu¬ 
lation (see (jUl) 

I F"'{U)Uldx = ^ F"'{U)UldU 

= 2 (f F'"{U)f{U)dU = 0 


for regular functions F(U) in the region of varia¬ 
tion of U. Thus, we can use unperturbed Whitham 
equations for Eq. (1551) and transform the results to 
Eq. (H51) by means of the substitution (|5i)). 

Let us illustrate this approach by its application 
to the Gurevich-Pitaevskii problem of evolution of 
initial step-like pulse 


u(x,0 ) 


u_ for X < 0, 
0 for X > 0, 


(56) 


according to Eq. (H51) . The substitution (15^ trans¬ 
forms this problem to the same problem for the 
equation (1551) with U- = U- — eF(u-). The well- 
known solution of the Gurevich-Pitaevskii problem 
for unperturbed Whitham equations applicable to 
(1551) yields, in particular, the speeds of edges of the 
dispersive shock wave (see, e.g., iSH), 


s_ = —6U- = —6(u- — €F(u-)), 
s+ = 4C7_ = 4('u_ — eF(u-)), 


(57) 


where s_ is the speed of the small-amplitude edge 
and s+ is the speed of the soliton edge. For 
F(u) = these formulae coincide with the results 
of exact theory developed in [s^ for the Gardner 
equation and it is instructive to compare them with 
the results of El’s method 11| applied to Eq. (HHl) . 
According to this method, the speeds of the shock 
edges are equal to (see section IV in Q) 




= 54!^, (58) 


where uJo{u,k) and ujs{u,k) are the linear and the 
‘soliton’ dispersion laws, respectively, given in our 
case by 

uJo(u, k) = V{u)k — k^, 

^ ~ ~ ~ ’ (59) 

ujsiu, k) = V{u)k -\- k^. 

and V{u) = 6{u — eF(u)). In the problem (1551) we 
have M+ = 0 and fc_ and kg should be found by 


solving the differential equations 


dk 

dcoo/du 

k{u+) = 0, 

du 

V{u) — dujo/dk' 

dk 

duJs/du 

1 

II 

o 

du 

V{u) — dujg/dk’ 


Easy integration yields 

k- = ks = 2\/u- — eF{u-) 


(60) 


(61) 


and substitution of these values into dSHl) gives 
Eqs. (|57)) . 

It is remarkable that the perturbation theory re¬ 
produces in this case the results correct beyond ap¬ 
plicability of the perturbation approach. Besides 
that, in framework of the perturbation theory under 
consideration the dissipative effects can be easily 
taken into account by adding corresponding terms 
to Eq. (|55)) and by calculation of their contribution 
into the right-hand sides of the Whitham equations 
(HI- One may suppose that this approach can be 
useful for consideration of the problems of the type 
considered in a recent preprint [34| where a com¬ 
bined action of nonlinear, dispersive and dissipation 
effects should be taken into account. 

At last, the perturbation theory is not limited 
to the step-like initial conditions. For example, if 
one wishes to consider evolution after wave breaking 
described by the Gardner equation (l43l) , so that the 
initial condition can be reduced to the form (see, 
e.g., [HI) 

X — 6{u — €u^)t = —u^ , (62) 

then the substitution (H51) transforms the problem 
to solving the equation (1551) with the initial condi¬ 
tion 


x-6Ut = -U'^ - 3et7^ 


(63) 


Here the Whitham equations correspond to the un¬ 
perturbed KdV case. Their solution with the initial 
condition (1551) for e = 0 was found by Potemin , 
and its generalization to the case e 0 can be done 
by including higher commuting flows of the KdV 
hierarchy into the standard method of integration 
of unperturbed Whitham equations (see, e.g., nil). 


3. Steady state solutions of perturbed 
Whitham equations 

In dissipative-dispersive systems one can distin¬ 
guish several characteristic stages of evolution of 
a pulse. For example, if we consider evolution of 
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a step-like initial pulse (I56L then it is natural to 
expect that the dispersion term with third order 
x-derivative is much more important at the initial 
stage of evolution than the dissipative perturba¬ 
tion term typically proportional to the second or¬ 
der derivative Uxx (or even u itself; see, e.g., (^ 1 . 
Hence, if we introduce small parameter R ^ e, 
then for time t < 1 /e we can neglect the dissipa¬ 
tive perturbation and evolution is described by the 
classical solution of the Gurevich-Pitaevskii prob¬ 
lem Q. However, for t > 1/e the damping of soli- 
tons in the dispersive shock wave train becomes es¬ 
sential and this damping can be compensated by 
the non-zero boundary condition at the small am¬ 
plitude limit X —>■ —00 where u —>■ U-. Thus, 
for t 1 /e we arrive at the steady-state solu¬ 
tion of the perturbed KdV equation with dissipa¬ 
tion balanced by the non-zero boundary condition. 
If damping is small enough, then the steady-state 
dispersive shock can be described by a modulated 
cnoidal wave solution (1261) where \i are the func¬ 
tions of ^ = X — Vt and satisfy the perturbed 
Whitham equations (HI- They are greatly sim¬ 
plified by noticing that Si = Ai -I- A 2 -I- A 3 is their 
integral si = const provided V = —2si and, hence, 
the Whitham equations can be reduced to 


dAj 


((2Ai - Si - u)R) 
8 ~ ^j) 


i = 1,2,3. 


(64) 


It is easy to check that these equations indeed have 
the integral si = const so that above ansatz is 
justified. Easy calculation shows that two other 
symmetric functions S 2 = A 1 A 2 -|- A 1 A 3 -I- A 2 A 3 and 
S 3 = A 1 A 2 A 3 (see dSni) satisfy the equations 


^ = ^ = + (65) 

Thus, we have reduced the problem to the system 
of two ordinary differential equations for S 2 , S 3 and 
Ai are considered as functions of S 2 , S 3 being the 
roots of the algebraic equation 


Instead of introduction of new dependent variable 
S 3 we can return in this case to Eqs. dH and 
consider, say, Ai and A 2 as functions of A 3 where 
A 3 = Az{^). Then we get 

dXi _ A3 — A2 dXi _ A3 — Ai . . 

^ “ A2 - Ai ’ ^ ““A2-A1’ ^ ^ 

and this system has, as we know, two integrals 


Ai -(- A 2 A 3 — Si — const, 

A 1 A 2 -l- A 1 A 3 -|- A 2 A 3 = S 2 = const. 


(69) 


Consequently, Ai and A 2 as functions of A 3 are two 
roots of the quadratic equation 

A^ — (si — A 3 )A -|- S 2 — (si — A 3 )A 3 = 0. (70) 

They are ordered according to inequality Ai < A 2 
and the integration constants si, S 2 can be found 
from the boundary conditions. Substitution of ex¬ 
pressions for Ai = Ai(A 3 ) and A 2 = A 2 (A 3 ) into dif¬ 
ferential equation for A 3 (^) (see (IMll l with known 
(uR) solves in principle the problem. However, 
some important consequences can be obtained with¬ 
out integration of this differential equation. Let us 
illustrate this by application of the derived formulae 
to the steady-state dispersive shock wave evolved 
from the initial step-like distribution (I56|) . 

The initial distribution (IMl) suggests that at x —>■ 
—00 the shock has a form of a small-amplitude wave 
with m —0, A 2 —>■ A 3 , so that here according to 
(1^ we have Ai = A/ = —U- and the integrals 
take the form 


si — — 'll— -(- 2 A 3 , 52 — — 2u—X^ -\- (Ag )^. (71) 

At the soliton edge m —>■ I, A 2 —)> Ai the wave (1^ 
reduces to 

^ ^ cosh2[VA3-Ai (x - Vt)] 

and since from (15611 we have the boundary condition 
u(x, t) — >■ 0 as X — >■ 00 , we get here 

A+=A+ A+=0. (72) 


A^ — siA^ -|- S2A — S3 = 0. (66) 

Especially simple and practically important sit¬ 
uation realizes when {R) = 0, hence we have the 
second integral S 2 = const and the problem reduces 
to a single differential equation 


Then the leading soliton at the soliton edge is de¬ 
scribed by the equation 

u(x,t) = - 5 — -. (73) 

cosh^[-v/—Ai (x -I- 4Ait)] 

Now we substitute the boundary conditions dUl) 
into Eq. dini) to obtain the relation 

sf - 4 s2 = 0 
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(74) 









between the integration constants which must be 
fulfilled along the whole dispersive shock wave. Its 
application to Eqs. (El) yields at the small ampli¬ 
tude edge A]" = —u-, A^ = A;^ = —it_/4 and 
hence 

= S2 = —u_. (75) 

Then equation (1701) gives at the soliton edge where 
Ag = 0 the value of a double root A)^ = A^ = 
—3u_/4 and as a result we find the expressions for 
the soliton amplitude Og = —2Ai and the velocity 
V = —2si of the shock wave in terms of a given 
value of the initial step amplitude, 

3 

Qs = 2 ^-’ ^ “ 3u_. (76) 


in any parametrization with account of perturba¬ 
tion terms and after that transform these equa¬ 
tions to the known Riemann invariants of the un¬ 
perturbed system arriving at perturbed Whitham 
equations. Thus, one may hope that the direct 
Whitham method can find in future many inter¬ 
esting and important applications. 
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